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Uncovering the mechanisms underlying robust re- 
sponses of cells to stress is crucial for our under- 
standing of cellular physiology. Indeed, vast 
amounts of data have been collected on transcriptional 
responses in Saccharomyces cerevisiae. However, only 
a handful of pioneering studies describe the dynamics 
of proteins in response to external stimuli, despite the 
fact that regulation of protein levels and localization is 
an essential part of such responses. Here we charac- 
terized unprecedented proteome plasticity by systemati- 
cally tracking the localization and abundance of 5,330 



yeast proteins at single-cell resolution under three differ- 
ent stress conditions (DTT, H2O2, and nitrogen starva- 
tion) using the GFP-tagged yeast library. We uncovered 
a unique "fingerprint" of changes for each stress and 
elucidated a new response arsenal for adapting to radi- 
cal environments. These include bet-hedging strategies, 
organelle rearrangement, and redistribution of pro- 
tein localizations. All data are available for download 
through our online database, LOQATE (localization and 
quantitation atlas of yeast proteome). 



Introduction 

The budding yeast adapts its ptiysiology to radical external 
and internal fluctuations in a rapid and robust manner. It sur- 
vives exposure to toxic substances, consumes a variety of 
nitrogen and carbon sources, lives both as a haploid or diploid 
organism, and carries out both mitosis and meiosis. Uncov- 
ering the mechanisms underlying its robustness in the face 
of extreme environments, genetic mutations, and life phases 
motivates many explorations of this eukaryotic cell. Indeed, 
much effort has been invested in measuring the transcriptional 
responses of whole yeast populations under varied conditions 
over time (Gasch et al., 2000; Hughes et al., 2000; Causton 
et al., 2001; Nagalakshmi et al., 2008). However, it is now 
widely appreciated that posttranscriptional, translational, and 
posttranslational regulation are also required for cell adapta- 
tion to stimuli (Stathopoulos-Gerontides et al., 1999; Beyer et al., 
2004; Hedbacker et al., 2004; Newman et al., 2006; Sigal 
et al. 2006, 2007; Erjavec et al., 2007; Aragon et al., 2009; 
Frenkel-Morgenstern et al., 2010; Eden et al., 2011; Lee 
et al., 2011; Tkach et al., 2012). Because proteins are the 
major components that carry out actions within the cell, and 
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are extremely sensitive to their surrounding milieu, even mild 
changes in protein abundance or subcellular localization may 
lead to drastic consequences in cellular physiology. We there- 
fore wished to add to previous efforts in studying physiologi- 
cal robustness by characterizing yeast stress responses at the 
protein level, focusing on two key properties: protein localiza- 
tion and abundance, at a single-cell resolution. 

Results 

A novel system for measuring 
genome-wide protein abundance and 
localization at single-cell resolution 

We have constructed a high-throughput microscopy setup to 
measure the levels of 5,330 yeast proteins and manually assign 
their localization using strains from the GFP fusion library (Huh 
et al., 2003). In this collection of strains, a GFP was integrated 
C-terminally into the genomic locus of each yeast ORF, main- 
taining the chromosomal context and preserving the natural 
promoter of each gene (Huh et al., 2003; Newman et al., 2006). 
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Figure 1 . Single-cell quantification of GFP-tagged protein levels reveals diverse proteomic changes during nitrogen starvation, (a) To track single cells, we 
integrated a cassette for expressing cytosolic mCherry into every strain of the yeast GFP collection, (b) We then automatically acquired images of hundreds 
of individual cells from each strain during growth in control conditions (logarithmic growth in SD) and three different stress stimuli. Using the mCherry (1), 
we segmented images (2) into single cell objects (3) and extracted GFP intensity in a single-cell resolution (4]. Manual inspection of the library enabled the 
determination of subcellular localization for each protein (bottom). Bars, 5 pm. (c) A scatter plot (f^ = 0.97] of median fluorescence values measured for 
each of the 5,330 strains of the yeast GFP library in two independent experiments in SD (a.u., arbitrary units), (d) Bubble plot depicting relative abundance 
(in SD) indicated by the size of the circle and sorted by the log2 ratio of all measured strains under nitrogen starvation relative to SD. Red represents up- 
regulated events, blue represents down-regulated events, and block circles mark strains that were found to be significantly nonunimodol (i.e., to exhibit 
two subpopulations) by a Hartigon's dip test, (e) Single-cell resolution analysis reveals bimodol distribution of expression of ribosome particles and ribosome- 
associoted chaperones during nitrogen starvation relative to SD. Insets shows mCherry expression to remain unimodal in each strain under the same condi- 
tion. This experiment was performed once on at least 1 GO cells per sample. 



To facilitate systematic quantification of GFP fluorescence in 
these strains, we introduced a constitutively expressed cytosolic 
marker into the entire library (77?F2pr-mCherry; Fig. 1 a; Tong 
et al., 2001; Cohen and Schuldiner, 2011). The cytosolic fluor- 
escence allows us to segment the images into hundreds of indi- 
vidual yeast cells for each strain and to extract protein abundance 
(detected as GFP fluorescence) under the desired condition 
(Fig. 1 b and Fig. SI). 

We first tested the reproducibility of our automated mi- 
croscopy approach for protein abundance. During logarithmic 
growth in synthetic medium with dextrose (SD), our measure- 
ments were highly reproducible (r^ = 0.97, slope = 0.98, P < 
0.01) in two biologically independent acquisitions of all strains 
(Fig. 1 c and Table S 1). To assay the accuracy of our system, we 
compared our data to previously published protein measure- 
ments of yeast grown in the same medium. We found high agree- 
ment (r^ = 0.79 and r = Q.ll) between our system and the levels 
of the 2,153 strains from the same GFP library previously mea- 
sured by flow cytometry (Newman et al., 2006) and measure- 
ments taken for the same library with an independent microscopic 



setup (Andrews, B., personal communication), respectively 
(Fig. S2). Our results also correlate with the levels of native, 
untagged, proteins measured by mass spectrometry (Walther 
et al., 2010), Western blotting on tandem affinity purification 
(TAP)-tagged strains (Ghaemmaghami et al., 2003), and ribo- 
somal footprint values (Ingolia et al., 2009; Fig. S2). 

In addition to measuring protein abundance, images were 
used to extract localization information by manual inspection 
(Fig. 1 b, bottom). We categorized each protein into one of 13 
subcellular localizations that can be distinguished without the 
use of colocalization markers. The main advantage of manual 
inspection is the miniscule error rate compared with faster auto- 
mated methods (Rimon and Schuldiner, 201 1). To enable accu- 
rate determination of a low-abundance cytosolic signal as bona 
fide expression (in contrast to cellular autofluorescence), we set 
the threshold of autofluorescence in each condition by measuring 
a library strain that does not harbor GFP (Fig. SI). Therefore, 
we could accurately assign expression and novel localizations 
to 72 proteins that are expressed at low abundance and were not 
visualized by previous, less sensitive microscopic platforms in SD. 
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We also identified 48 proteins expressed in a condition-specific 
manner. 21 of the newly visualized proteins are essential, and so 
reveal new important information on proteins that are vital for 
yeast cellular life (Fig. SI and Table S2). 

Intensity analysis reveals trends 
of abundance cl-ianges under stress 

After establishing the accuracy and reproducibility of this meth- 
odology, we used our platform to track proteomic changes dur- 
ing growth in three stress conditions that are known to activate 
distinct cellular responses: oxidative stress caused by hydrogen 
peroxide (H2O2), reducing stress caused by DTT, and nitrogen 
starvation (Gasch et al., 2000; Travers et al., 2000; Giaever et al., 
2002; Hillenmeyer et al., 2008; Tyedmers et al, 2008). We found 
that at a single measured time point each condition evoked hun- 
dreds of changes in protein abundance (Fig. S3; for the full lists 
of up/down-regulated proteins, see Tables S3, S4, and S5). These 
changes were indicative of the responses required to restore ho- 
meostasis in this condition. For example, during nitrogen star- 
vation, 336 proteins were present at a significantly higher level 
(P < 0.01; Fig. 1 d). These proteins were enriched for involvement 
in oxidation-reduction processes (P = 2.8 x 10"'') and hetero- 
cycle metabolic processes (P = 2.9 x 10~'). In addition, the 360 
proteins that were present at significantly lower levels (P < 0.01) 
were enriched for ribosomal RNA (rRNA) processing (P = 2 x 
lO"''') and nitrogen compound metabolic process (P = 1.4 x 10"^^). 
Importantly, comparing these direct proteomic measurements 
to microarray data revealed that over a third of the observed 
changes in protein levels (in all three stresses) could not be pre- 
dicted by changes at the mRNA level (Fig. S4 and Tables S3, S4, 
and S5; Greenbaum et al., 2003; Newman et al., 2006; Walther 
et al., 2010; Schwanhausser et al., 2011). This suggests a much 
more dramatic, proteome-wide, posttranscriptional regulation of 
protein abundance levels than previously reported. 

Measuring protein abundance using single-cell resolution 
allowed us to explore variations in protein abundance under stress 
beyond the population level described above and beyond that 
possible by the pooled samples required for microarray analysis. 
For example, we detected 57 proteins that acquired a bimodal 
distribution upon nitrogen starvation (P < 0.05; marked as black 
circles in Figs. 1 d and 6; representative examples are shown in 
Fig. 1 e; for the full list of proteins that changed modality under 
stress see Table S6). This response is not a reflection of cell size 
(not depicted) or genotype because the mCherry expression in 
these cells remained unimodal (Fig. 1 e, insets). The bimodal re- 
sponsive proteins were enriched for subunits of the large ribo- 
somal particle (P < 7 x 10~^) and cellular protein metabolic 
processes (P < 0.001), which suggests that populations of yeast 
cells distribute their responses such as to reach two opposite out- 
comes: cell proliferation versus survival (Joo et al., 201 1). 

Single cell analysis uncovers a bet-hedging 
strategy in response to stress 

Recent studies observed phenotypic heterogeneity in altered 
environmental situations displayed by isogenic bacterial and 
yeast populations (Becskei et al., 2005; Allen et al., 2006; Acar 
et al., 2008; Gefen et al, 2008; Benbadis et al., 2009; Gefen 
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Figure 2. A bet-hedging strategy in yeast underlies prolonged survival 
under starvation, (a) Flow cytometry dot plots show the transition of uni- 
modal distributed population of Pre3-GFP in SD (left] into two distinct 
subpopulations after 24 h of starvation (right; a.u., arbitrary units). High 
and low fractions (marked with broken elNpses on the right plot) were 
sorted accordingly. Microscopic analysis of the two sorted subpopulations 
demonstrates distinct morphologies. Bars, 5 pm. (b) Growth of the two 
subpopulations in SD immediately after sorting (24 h starvation, left) versus 
4 d after sorting (120 h starvation, right) of both Pre3-GFP and Ssbl-GFP 
(OD). Each growth well contained 100,000 cells at the starting point. 
These measurements are the representative means of two repeats; error 
bars indicate standard deviation in black. 

and Balaban, 2009; Halfmann et al., 2010; Meyerovich et al., 
2010; Zhang and Rainey, 2010; Balaban, 2011; Levy and Siegal, 
2012; Levy et al., 2012). This phenomenon, termed bet hedg- 
ing, specifically refers to survival strategies where each sub- 
population performs differentially depending on the type of 
environment that becomes selective. For a population, this pro- 
vides a long-term fitness advantage (de Jong et al., 201 1). Because 
our observation of bimodality could underlie a "risk-spreading" 
strategy, either attenuating or augmenting the cell's ability to 
survive under starvation, we tested the hypothesis of bet hedg- 
ing for two representative proteins: Pre3 (the (31 subunit of 
the 20S proteasome) and Ssbl (a ribosome-associated molecu- 
lar chaperone). In strains expressing the GFP-tagged proteins, 
we sorted out the low and high GFP-expressing fractions 
after 24 h of nitrogen starvation (Fig. 2 a, middle). By light 
microscopy we could observe that the cells in each fraction dis- 
played unique morphological characteristics: the high GFP- 
expressing fraction is similar to logarithmically growing cells, 
whereas the low GFP-expressing fraction is granular and re- 
sembles the quiescent fraction in stationary phase culture (Allen 
et al., 2006; Benbadis et al., 2009; Fig. 2 a). We then measured 
the growth rate of the two subpopulations (100,000 cells in 
each) either after immediate transfer to nutrient-rich conditions 
(SD) or after prolonged starvation (an additional 4 d) followed 
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by return to nutrient-rich medium. Indeed, we found that the 
high-GFP fraction had a dramatic advantage over the lower one 
(high-GFP fractions grew similarly to nonstarved culture; not 
depicted) if return to rich medium was rapid. However, the low- 
GFP fraction had clear fitness advantage over the high-GFP 
fractions if starvation was prolonged (Fig. 2 b). 

Integrating localization with abundance 
data reveals stress-specific whole 
organelle rearrangements 

Spatial organization in the organelle-bearing eukaryotic cell is 
highly regulated and important for cell viability. A fundamental 
aspect of stress responses is therefore how protein abundance 
changes in the context of spatial organization. To track the ex- 
tent of changes in the protein repertoire of each organelle, we 
calculated the mean log2 fold change of all proteins that were 
localized to every organelle under each condition (Fig. 3 a). 
Using this strategy, we observed that each stress condition pref- 
erentially affected specific organelles, either increasing or re- 
ducing their protein content. 

For example, under nitrogen starvation, where protein 
translation is reduced, we saw dramatic shrinkage of the nu- 
cleolus, most likely due to a reduction in rRNA production 
(Fig. 3, a and b). Under the same condition, plasma membrane 
content rose, most probably as a result of the cells' attempt to in- 
crease the presence of transporters on the cell surface and thus 
enhance uptake of sparse nutrients. An example of such an 
up-regulation is the family of allantoin transporters, which are 
required for uptake of alternative nitrogen sources (Chisholm 
et al., 1987; Isnard et al., 1996; Kaur and Bachhawat, 2007), 
that were below detection threshold under control conditions 
but became highly expressed when cells were starved for ni- 
trogen (Fig. 3 c). As another example, after DTT treatment we 
observed up-regulation of actin filament subunits (P = 0.002). 
This is intriguing in light of previous publications that sug- 
gested a major role for cytoskeletal filaments in induction of 
the secretory pathway stress response: the unfolded protein 
response (Aragon et al., 2009). 

We chose to focus on the changes that occur under H2O2 
treatment, where we observed a prominent deflation of the mito- 
chondrial proteome. This was accompanied microscopically 
by fragmentation of mitochondria (Fig. 3 e). To determine 
whether these changes are required to combat H202-induced 
stress, we first deleted proteins responsible for mitochondrial 
fission (Mozdy et al., 2000) and verified that in their absence 
mitochondria cannot fragment even under oxidative stress. These 
deletion strains indeed show a growth defect in the presence of 
H2O2 (Fig. 3 f). To broaden our perspective of essential func- 
tions of mitochondrial proteins while combating H202-induced 
stress, we focused on several nonessential mitochondrial proteins 
that were found to be essential for growth in H2O2 (Hillenmeyer 
et al., 2008). Following the fluorescent tagged fusions of these 
proteins by time-lapse resolution, we could find examples of 
proteins whose fluorescence signal changed dramatically as 
early as 30 min from initiation of the stress. Overlay with bright- 
field images suggests that these changes represent presence on 
the vacuolar membrane (Afg3, Sdh2, Mss51, Cbp3, Mdl2, and 



Aim37) or in the vacuole lumen (Mrplll, Pkpl, Sws2, Mrpl23, 
and Ymcl; Fig. 3 d). Interestingly, this phenomenon occurred for 
only '^30% of cells in the population. However, in those cells 
where it was seen, it was accompanied by accumulation of a 
matrix-targeted dsRed (MTS-dsRed) in the vacuolar lumen (as 
can be seen by overlay with bright-field images; Fig. 3 d) and by 
degradation of the protein over time (Video 1). Important to note 
is that most mitochondrial proteins (341/353) do not change 
their subcellular location during oxidative stress, providing a 
robust and inherent control for overall mitochondrial morphology 
under these conditions (Video 2). However, it remains formally 
possible that under oxidative stress mitochondria may swell up in 
a fashion that is so drastic that they may resemble vacuoles in 
bright-field images. 

Three of the proteins whose localization were altered 
(Sws2, Mrpl23, and Mrplll) are components of mitochondrial 
ribosomes, therefore their elimination would enable a major and 
rapid down-regulation in mitochondrial proteome production. 
This sheds new light on how cells clear dysfunctional mitochon- 
dria in a timely manner under oxidative stress and ensure rapid 
reduction of mitochondrial activity. The full list of proteins dis- 
playing an altered pattern of their predominant localization under 
H2O2 is enriched for proteins essential under this condition (35%, 
P = 0.001; Fig. S5), which suggests that localization change is an 
integral component of the response arsenal to combat oxidative 
stress (for the responsive proteins that are essential under starva- 
tion but not in reference growth condition, see Fig. S5). 

Repatterning of subcellular localization 
occurs during exposure to stress 

Previous literature demonstrated the importance of localization 
changes as a regulatory strategy in a handful of specific proteins 
(Cyert, 2003; Erjavec et al., 2007). Using our platform, we set 
out to explore these phenomena on a proteomic scale. In addi- 
tion to expression levels, we systematically determined protein 
localizations under the same three stress conditions relative to 
the reference growth condition. Unexpectedly, a large number 
of proteins (235) displayed a change in their predominant sub- 
cellular localization under stress (Figs. 3, 4, and S5). Changes 
included nuclear proteins that could now be visualized in the 
cytoplasm and vice versa, plasma membrane proteins that could 
now be visualized in the vacuole, pan-ER localization that was 
visualized as nuclear vacuolar junctions, mitochondrial proteins 
that could be visualized in the vacuole, and proteins from all 
organelles that could be visualized in punctate foci (for the full 
scheme of changes observed see Fig. 4 d and for full list of 
changes see Table S7). Indicative of the functionality of the 
changes in protein localizations is that at a very high occurrence 
(48%) we found concerted changes for physically interacting 
proteins (interacting proteins were determined by Collins et al., 
2007; for full details see Table S8). The unprecedented extent 
and diversity of changes demonstrate the unappreciated capac- 
ity of a wide variety of proteins to shift localization. 

Because protein-protein interactions mediate cellular 
activity, we hypothesized that the changes that we observed 
in protein localization may be accompanied by a change in 
the networks of physical interactions. To test this, we mapped 
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Figure 3. Integration of proteome abundance and subcellular organization enables characterization of the dynamics of organelle morphology and com- 
position under stress, (a) Hierarchical clustering of mean log2 ratio of the change in abundance of proteins in each organelle under stress relative to SD. 
(b] Representative images of nucleolus proteins in SD or during growth in starvation medium. Colocalization with a nuclear mCherry (NLS) allows visualiza- 
tion of the massive shrinkage of this compartment relative to the nucleus under these conditions, (c) Representative images of cell periphery proteins that are 
expressed only under starvation and cause an increase in the protein levels in this organelle, (d) Representative images of mitochondrial proteins tagged 
with GFP coexpressed with a matrix-targeted dsRed (MTS-dsRed) during growth in SD or hl202 treatment. During stress these proteins become localized to 
the vacuole lumen or membrane (as can be seen by the corresponding bright-field image) concomitantly with presence of the mitochondrial dsRed signal 
in the vacuole, (e) Mitochondrial fragmentation occurs when cells are exposed to medium containing 1 mM H2O2. This fragmentation is blocked in the 
absence of the fission machinery proteins Fisl , Dnm 1 , or Mdvl . (f) Fragmentation seems to be required for growth during oxidative stress as loss of fission 
capacity reduces fitness of cells under 1 mM H2O2 relative to control (YEPD). Bars, 5 \im. 



the entire physical interaction network for two proteins that 
changed their localization under DTT stress using the protein 
fragment complementation assay (Tarassov et al., 2008). The 
first, Pin3, is a protein of unknown function whose overproduc- 
tion induces the appearance of prions (Derkatch et al., 2001). 



Pin3-GFP shifted from homogenous cytosolic distribution to 
discrete punctate foci. Indeed, under DTT Pin3 acquired four 
new significant interactions with the proteins Arcl5, Arcl8, 
Lasl7, and Thsl (P < 10"'^). Arcl5, ArclS, and LaslV are all 
part of cortical actin patches that would give a punctate foci 
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Figure 4. Protein localization is dynamic under stress, (a) Two strains per condition (out of a total of 254 unique changes in protein localization) as 
examples for changes in predominant protein localization between growth in SD and growth under stress, (b] The only localization change that is common 
under all three stresses is proliferation of P-bodies. Shown ore two P-body proteins, Dcpl and Dcp2, changing localization from the cytosol to punctate 
foci. Broken lines represent the contour of the cell from bright-field images, (c) Foci of Dcpl and Dcp2 colocalize with Edc3-mCherry during growth under 
every stress condition, which shows that these are bona fide P-bodies. (d] A schematic diagram depicting the types and numbers of changes in cellular 
localization observed during yeast growth in the three different stresses. Bars, 5 pm. 
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Figure 5. Examples of relocalizations of proteins from the cytosol into the nucleus during growth In nitrogen starvation medium. 71 cytosolic proteins 
categorized to eight different cellular functions could be visualized in the nucleus during starvation, (a) A representative protein from each of the categories 
is shown in SD and under nitrogen starvation, colocalized with a nuclear (NLS)-mCherry. Bars, 5 fjm. (b) Pie chart with the relative size of each functional 
group and information on the genes categorized in it. 



pattern by fluorescent microscopy. In support of these new 
interactions, the shift of Pin3 to interact with the actin cyto- 
skeleton under stress has been recently demonstrated (Madania 
et al., 1999; Ganusova et al., 2006; Chernova et al., 2011), 
validating our hypothesis that changes in protein localization 
may reflect the acquisition of new cellular functions. The sec- 
ond, Accl, catalyzes the carboxylation of acetyl-CoA to form 
malonyl-CoA (Hasslacher et al., 1993). Interestingly, during 
DTT treatment its predominant localization changed from the 
cytosol to rod-like foci. In accordance with this movement, 
Acc 1 acquired novel interactors after DTT treatment (for a full 
list see Table S9). Of note are the new physical interactions 
with both Pre9 (the a3 subunit of the 20S proteasome) and 
Rpn8 (the regulatory subunit of the 26S proteasome), which 
suggests the involvement of Accl in proteasomal degradation 
during exposure to DTT. 

However, by far the most remarkable change in localiza- 
tion during starvation was observed for 7 1 cytosolic proteins 



whose predominant subcellular distribution became nuclear 
(Fig. 5). Although six of these proteins are transcription fac- 
tors, whose change can be easily rationalized, the vast major- 
ity are proteins not previously known to be required in the 
nucleus under such conditions. For example, eight proteins 
are chaperones and eight are oxido/reductive homeostasis- 
related proteins. A striking observation was that 20 meta- 
bolic enzymes entered the nucleus (Fig. 5). It has been shown 
recently that the two mitochondrial enzymes fumarase and 
aconitase play a minor function in the nucleus that becomes 
essential under DNA-damaging conditions (Yogev et al., 
2010; Ben-Menachem et al., 2011). Our findings would indi- 
cate that such cases may be much more prevalent than pre- 
viously appreciated. 

Most localization changes were condition-specific (Table S7). 
However, four proteins changed localization in a concerted 
manner during all three stresses (Dcpl, Dcp2, Keml, and Edel ; 
Fig. 4 b). Three out of the four proteins are localized to P-bodies 
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Figure 6. Infegration of all proteomic changes gives a "thumbprint" of 
cellular stress responses. Representative schemes integrating all proteomic 
changes per each stress (a, starvation; b, H2O2; c, DTT) demonstrate the 
unique thumbprint of cell remodeling in each of the responses. The entire 



(Fig. 4 c) and play essential functions in P-body formation, sup- 
porting the previously published observation that P-body func- 
tion and proliferation play a major role during diverse stresses 
(Nissan et al., 2010). 

Collectively, our data open a new porthole to the vast 
extent of changes exhibited by the cell under stress solely at 
the level of the proteome. A birds-eye view of these changes 
can be seen when integrating localization and abundance data 
into unique "fingerprints" of the yeast cellular response to each 
stress (Fig. 6). 

□iscussion 

The remodeling of expression patterns that occurs while adapt- 
ing to new environments was characterized extensively in the 
past decade using first microarray and now deep-sequencing 
technologies. This fundamental information has provided sig- 
nificant clues about the way in which we understand stress 
responses today. Interestingly, previous publications demon- 
strated only a small overlap between the groups of genes es- 
sential to combat a specific stress and the mRNA transcripts 
that change under the same conditions (Giaever et al., 2002; 
Hillenmeyer et al, 2008). This strongly suggests that proteomic 
level changes in stressful environments may be predominant 
and can give us new insights about the cellular state. Indeed, 
accumulated literature of proteomic studies promises that there 
are many more mysteries to be elucidated regarding the cellu- 
lar response to environmental stress. Using a high-throughput 
microscopy system, we show here the first nonbiased system- 
atic high-resolution analysis of protein abundance and local- 
ization changes at single-cell resolution after three different 
stress conditions. Moreover, we integrate these two types of 
data to give a fundamental perspective of abundance changes 
within the context of spatial eukaryotic organization to detect 
whole organelle proteomic composition and morphology for the 
first time. We have observed these stress-specific organelle 
changes to be crucial for cellular robustness in the face of stress. 

One notable example for how such data can be used is in 
the modeling of transcriptional networks that govern genome 
expression changes. Our work shows that there are extensive 
changes in the predominant localization of transcription factors 
under stress. This strengthens the notion that in many cases 
transcription factors are controlled not by their levels of expres- 
sion but rather by their localization. This would suggest that 
previous analyses that have used only mRNA levels to try and 
identify transcriptional networks are missing an essential aspect 
of transcription factor regulation. 

There are numerous advantages to probing cellular pro- 
cesses at a single-cell resolution. As an example, it allowed 
us to uncover the proteins that are regulated in a bi-modal 



proteome is divided into organelles represented as colored segments on 
the circumference. The black bar plot above each organelle gives each 
protein's abundance under SD and is sorted in each organelle. Blue (down- 
regulation] and red (up-regulation) bars demonstrate protein abundance 
changes. Lines between organelles demonstrate localization changes and 
are colored according to the destination organelle. 
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fashion in response to stress. Bet hedging is poorly under- 
stood at the molecular level, though it plays a major role in 
diseases, such as cancer (Levy and Siegal, 2012) and per- 
sistence of microbial infections (Levy and Siegal, 2012). 
Measuring the changes in protein distribution within the 
population under a variety of conditions should shed novel 
light on the molecular mechanisms used by species to adapt 
toward radical environments and promises future solutions to 
the current failure of chemotherapies and antibiotics on bet- 
hedging populations. 

The platform we present here may be used to measure 
the proteome under a variety of stress conditions, life phases, 
genetic backgrounds, and drug treatments. To spearhead these 
efforts and to facilitate easy access to our data we have cre- 
ated a database of changes in protein levels and localizations 
(LOQATE, localization and quantitation atlas of the yeast pro- 
teome) that will house the information from this current work 
as well as our future screens (www.weizmann.ac.il/molgen/ 
loqate). Such systems-wide, precise localization and quantita- 
tion, directly at the protein level and at single-cell resolution, 
opens up new avenues in the postgenomic era, which hold the 
promise to shed light on novel aspects of cell biology that have 
not been previously characterized. 

Materials and methods 

Strain management 

Insertion of cylosolic mCherry to the GFP library. Synthetic genetic array (SGA] 
technique was performed between a MATa haploid strain harboring TEF2p,- 
mCherry::URA3 integrated into the URA3 locus (plasmid for creation of the 
strain was a kind gift from D. Breslow, Stanford University, Stanford, CA; 
Breslow et al., 2008) against the GFP collection (::HIS3; the library was 
a kind gift from J. Weissmon, University of California, San Francisco, San 
Francisco, CA; Huh et al., 2003]. Mating was performed on rich media 
plates, and selection for diploid cells was performed on plates lacking 
both HIS and URA. Sporulation was then induced by transferring cells to 
nitrogen starvation plates for 5 d. Haploid cells containing all desired mu- 
tations were selected for by transferring cells to plates containing all selec- 
tion markers alongside the toxic amino acid derivatives Canovanine and 
Thiolysine (Sigmo-Aldrich) to select against remaining diploids and lacking 
Leucine to select for only spores with an "a" mating type (Tong et al., 2001 ; 
Cohen and Schuldiner, 201 1). SGA procedure was validated by inspect- 
ing representative strains for the presence of the GFP-tagged strains and for 
the cytosolic mCherry expression. To manipulate the collection in high-density 
format (384), we used a RoToR bench top colony arrayer (Singer Instru- 
ments). Yeast strains used in this study are in Table SIO. 

Strain grov^h. The manipulated 5,330 strains were grown in 50 ^jl 
SD (0.67% yeast nitrogen base without amino acids [Conda Pronadisa] 
and 2% dextrose) containing the appropriate supplements for selection in 
384-well plates (catalogue No. 781 162; Greiner Bio-One). 

Application of stress. For hydrogen peroxide (H2O2) treatment, cells 
were grown to early-log phase, at which point H2O2 (catalogue No. 
2186-01; J.T.Baker) was added at a final concentration of 1 mM. After 
60 min of incubation, plates were taken for microscopic imaging. 

For nitrogen starvation, cells were grown to early log phase and 
centrifuged mildly (1,000 g), then fresh SD (0.67% yeast nitrogen base 
without amino acids and without ammonium sulfate [Conda Pronadisa] 
and 2% dextrose) was added. After 15 h of incubation, plates were taken 
to microscopy imaging. 

For DTT treatment, cells grown overnight were back diluted into 
2 mM DTT containing SD medium. After 3 h of incubation, plates were 
taken for microscopic imaging. 

Plosmids and deletions. Deletion strains were prepared by replace- 
ment of the ORFs with a pCgMET15 cassette using homologous recom- 
bination with 40 bp of homology (Kitada et al., 1 995). pRS4 1 6 plasmid 
expressing MTS-dsRed under the ADHl promoter was provided by 



J. Nunnari (University of California, Davis, Davis, CA; Meeusen and 
Nunnari, 2003). pRS426 plasmid expressing NLS-tdTomato under the 
GPDl promoter was provided by D. Kaganovich (Hebrew University, 
Jerusalem, Israel; Kaganovich et al., 2008). Plosmids used in this study 
ore in Table SI 1 and primers used in this study are in Table SI 2. 

Automated imaging 

High-throughput fluorescence microscopy. Microscopic screening was per- 
formed using an automated microscopy setup as described previously 
(Cohen and Schuldiner, 201 1). Cells were moved from agar plates into 
liquid 384-well polystyrene growth plates using the RoToR arrayer. Liquid 
cultures were grown overnight in SD medium in a shaking incubator 
(LiCONiC Instruments) at SOX. A JANUS liquid handler (PerkinElmer), 
which is connected to the incubator, was used to bock-dilute the strains 
to ~0.25 CD into plates containing the same medium. Plates were then 
transferred back to the incubator and were allowed to grow for 3.5 h 
at 30°C to reach logarithmic growth phase, as was validated in prelimi- 
nary calibration. The liquid handler was then used to transfer strains into 
glass-bottom 384-well microscope plates (Motricol Bioscience) coated 
with Concanavalin A (Sigma-Aldrich) to allow cell adhesion. Wells were 
washed twice in medium to remove floating cells and reach cell mono- 
layer. Plates were then transferred into an automated inverted fluorescent 
microscopic ScanR system (Olympus) using a swap robotic arm (Hamil- 
ton). Imaging of plates was performed in 384-well format using a 60x 
air lens (NA 0.9) in SD medium at 24°C with a cooled charge-coupled 
device camera (ORCA-ER; Hamomatsu). Images were acquired at GFP 
(excitation at 490/20 nm, emission at 535/50 nm) and mCherry (excita- 
tion at 572/35 nm, emission at 632/60 nm) channels. 

Image analysis. Our screening assay was designed to explore yeast 
cell biology by assessing two cellular key features of interest: subcellu- 
lar localization and fluorescence intensity. To analyze these images we 
have used an in-house script to browse manually and assign localization 
rapidly and efficiently. To extract proteomic abundance from images, we 
used the Olympus ScanR analysis software. This allows for the preprocess- 
ing of images by background subtraction, and segmentation of images to 
identify individual cells as separate objects. Specifically, we performed 
the following steps: 

(a) Segmentation on the basis of the edge module of the cytosolic 
mCherry protein expression. 

(b) Background correction using the rolling ball algorithm. 

(c) Definition of measured populations. Because several measure- 
ments are collected for each cell (e.g., fluorescence intensity, area, shape), 
we have created a multiparameter gate to ensure that our population was 
homogenous and that data arise from clearly defined cells only. The mean 
GFP intensity for each object (cell) of each strain was extracted to Excel 
files (Microsoft), allowing data processing for single-cell resolution from 
within a given population. 

Data processing 

Median measurement. The median GFP intensity for each strain was mea- 
sured from the remaining objects for each strain under each condition after 
removing dead cells. Dead cells ore highly fluorescent and must be re- 
moved from the analysis because they artificially raise the mean GFP inten- 
sity values. Because the fluorescence and shape of dead cells fall within 
range of those features for normal cells, the software could not gate them 
out. To automate the removal of dead cells detected as objects by the 
ScanR software, for each strain we removed any objects that were very 
high outliers in their mean GFP measurement. Therefore, any objects with 
mean GFP measurement above UQ -1- 3 x IQR were removed from analysis, 
where UQ is the upper quartile and IQR is the interquartile range. Overall, 
>94% of cells originally screened were classified as alive. On average, 
5% of cells per strain were dead under stress. Because we performed two 
independent measurements under reference conditions (SD), we combined 
scores from both measurements to get one median and standard deviation 
value for each strain under reference conditions (summarized in Table SI). 

Removing strains from the analysis. The following strains were removed 
from the analysis: 

(a) Strains in which <25 objects were detected. 

(b) Strains whose proper subcellular localization has been demon- 
strated to be dependent on the C terminus of the protein. A full list of the 
systematically mislocalized proteins that were removed from analysis is 
given in Table SI 3. 

(c) Strains whose tagged ORFs are located near the CANl, LYPl, 
or URA3 loci. Such strains could not pass the SGA required to moke the 
Tef2-Cherry background. A full list is given in Table SI 3. 
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(d) Contaminated strains. Any strains showing localization different 
than that annotated in SGD or that shown in the original Huh et al. (2003) 
published dataset were assumed to be the result of contamination and 
were removed. 

(e) Bud neck proteins, as the analysis program could not accurately 
detect their signal, sometimes out of the detected cell boundaries. 

Detection of autofluorescence. Yeast cells not tagged with GFP emit 
fluorescence at some baseline intensity. We hove found this value to differ 
significantly across conditions. This may be a result of different cellular 
conditions under various stresses or of fluctuations in the light source. To 
account for this, we screened 86 randomly placed wells containing a wild- 
type strain with no GFP tag in each measured condition. After removing 
dead cells, the distribution of the intensity values for wild-type objects was 
not rejected for normality using the Shopiro-Wilk test (P < 0.01). We then 
approximated the autofluorescence as a normal distribution with mean 
autofluorescence value A„„difion ond standard deviation (Tco„diiion- Any 
strain whose median GFP intensity under a given condition was more than 
Acondifion + 2.58cr„„difi„„ is >99% likely to hove biological GFP expression. 
Strains below this cutoff are marked as "below threshold." 

If a strain falls below this threshold, but localization other than "cyto- 
sol" or "ambiguous" was assigned, it is kept in the analysis because we 
could be certain that we visualized that protein. 

Experimental controls. Testing the reproducibility. To test the stability 
of our microscopy platform, we plotted two independent measurements of 
all strains (5,330) in SD (Pearson's correlation coefficient on logarithmic 
scaled data tested = 0.97, slope = 0.98, P < 0.01). 

Measurement of accuracy. In order to test accuracy of our measure- 
ments, we plotted the calculated medians by our method versus two mea- 
surements of the same strains by flow cytometry and with the measurements 
made by a similar high-content screening setup, and found high reproduc- 
ibility (Fig. S2). 

Further testing. In order to show the relevance of these measurements to 
other proteomic datasets, we plotted the calculated medians by our method 
versus three measurements of yeast proteomes: native, untagged proteins mea- 
sured by mass spectrometry. Western blotting on TAP-tagged strains, and ribo- 
somal footprint values, and again found good agreement (see Fig. S2). 

Determining abundance change events 

For each condition, we determined which strains are significantly up- or 
down-regulated compared with their reference abundance levels. 

Preprocessing. In addition to the corrections described in "Data pro- 
cessing," we removed any strains that changed localization under the con- 
dition of interest compared with reference medium. Because the GFP signal 
is sensitive to cellular conditions (such as ionic strength or pH), levels can- 
not be compared between two proteins showing different localizations. 

Some proteins are only expressed under certain stress conditions 
and are not detected under reference conditions. To not miss these pro- 
teins, any strains that fell below the autofluorescence threshold in the refer- 
ence condition and therefore were not assigned localization, but were 
detected under the stress condition, were left in the analysis. Similarly, 
strains detected under reference but below the autofluorescence threshold 
under treatment were also included. 

Normalization of fluorescence signal. To make fluorescence values com- 
parable across reference and stress conditions, we normalized median fluor- 
escence values for each strain using the normalize. quantiles function in the 
R preprocessCore library (Bolstad et al., 2003). This method adjusts fluores- 
cence values such that the median values for each strain follow the same dis- 
tribution under each condition. The fluorescence values of objects for each 
strain under each condition are scaled to have the corrected median value. 

Definition of significant abundance change. All strains under SD in two 
independent experiments were sorted by their logio abundance values and 
binned such that each bin contained 5% of strains. For each bin, we plotted 
the distance from the diagonal for each strain in the scatterplot of replicate 1 
versus replicate 2. The distances within each bin were normally distributed 
(Shapiro-Wilk test, P < 0.01 ). For each stress condition, we plotted the logio 
abundance of each strain versus the overage of the logio abundance across 
the two reference conditions. We then determined the distance of each strain 
from the diagonal. Based on the distribution of distances between the two SD 
experiments in the corresponding abundance bin for each strain, we deter- 
mined an empirical P-value. Strains with P < 0.01 were marked as showing 
a significant abundance change (Fig. S3), and these events are marked 
above (red) and below (blue) the broken lines in Fig. S3. 

Fold change values. For each stress condition, we calculated the ratio 
of median GFP intensity measured under stress to the median GFP intensity 
under reference conditions. Ratios are given in Table SI . 



Testing multimodaliiy. In the case that not all cells are uniformly af- 
fected by a stress condition, the resulting distribution of fluorescence intensi- 
ties could be multimodal. The significance test described is not informative 
in this situation, and gives no indication that there might be multiple levels 
of response throughout a population of cells. To test each strain under 
each condition for multimodality, we used Hartigan's dip test (Hortigon 
and Hortigon, 1985), which tests how well the distribution of fluorescence 
for the objects fits the tightest fitting unimodal distribution. Strains in which 
unimodolity was rejected with P < 0.05 were marked as multimodal. This 
test was implemented using the dip. test command in the R dip test library. 
(The list of strains predicted to have multimodal distributions is given in 
Table S6). Manual inspection confirmed that all strains rejecting unimodol- 
ity were indeed bimodal. 

Comparison to published datasets 

Comparison to essentiality data. We used published essentiality datasets 
(Hillenmeyer et al., 2008) to determine which ORFs are essential under 
each condition tested. We downloaded the homozygous fitness defect 
scores from the web supplement to Hillenmeyer et al. (2008) and created a 
list of ORFs for each condition shown to be essential with P < 0.05. We used 
the conditions 5 mM H2O2 and SD for H2O2 and starvation, respectively. 

Comparison to mRNA data. For each ORF under each condition 
tested, we determined whether it was up-regulated, down-regulated, or 
showed no change at both the protein and mRNA level. Proteins were 
determined to be up or down-regulated based on criteria described in 
"Definition of significant abundance change." To determine if an ORF is 
up- or down-regulated at the mRNA level, we used published expression 
datasets. For DTT, we compared with the DTT time points at 1 5, 30, 60, 
and 120 min (Trovers et al., 2000). For starvation and H2O2, we com- 
pared with the nitrogen depletion and 0.32 mM H2O2 time points from 
30 min to 1 d and 1 0 min to 1 60 min, respectively (Gosch et oL, 2000). 
An ORF was defined as up-regulated if any time point showed a twofold 
increase in expression level compared with reference and as down- 
regulated if any time point showed a twofold decrease compared with 
reference. ORFs that hod time points showing both a twofold increase 
and twofold decrease were marked as ambiguous and were removed 
from the analysis. 

We then used these up/down/no change classifications to find 
ORFs marked as no change in mRNA levels but that did significantly 
change at the protein level according to our analysis. 

Protein-protein interactions analysis 

Dihydrofolate reductose-bosed protein fragment complementation assay 
(DHFR PCA). The assay for the yeast DHFR PCA followed entirely the pub- 
lished protocol (Tarossov et oL, 2008). In brief, A'lATa strains with the 
ORFs of Accl and Pin3 fused C-terminally to F[l,2] were mated to the 
entire MAJa collection of ORFs tagged with F[3]. The resulting diploids 
were subsequently selected for growth in the presence of methotrexate for 
positive DHFR PCA reconstitution with or without the addition of 3 mM 
DTT for 5 d in SOX. 

Data acquisition, colony quantification, and statistical analysis. Com- 
plete acquisition and analysis of each plate proceeded as follows. First, 
images of the diploid methotrexate selection with or without DTT were 
taken after 1 20 h of growth. Plate images were saved in JPG format at a 
resolution of 300 dpi. Using the freely available Bolony software (http:// 
code. google. com/p/bolony/downloods/list), the first step of the image 
analysis was to determine expected centers of the colonies arrayed in 
48 columns and 32 rows. To adjust for variation in plating and possible 
rotation of an image during image acquisition, we manually defined the co- 
ordinate center of a first colony in a first row of the array and of the lost 
colony of the last row. Results of positional array adjustment and detection 
of colony centers were also validated manually for all images. We then 
extracted the area for each colony and set a threshold of a positive inter- 
action to be >1 50. This threshold was chosen accordingly: in the negative 
interactions, distribution concentrated normally with a mean = 41, stan- 
dard deviation = 11.1, whereas positive interactions with controls distrib- 
uted with a mean = 400, standard deviation = 208. Therefore, each 
colony got a z-score according to its area and distance from the negative 
controls mean, and a P-value was calculated following the multiple hy- 
pothesis correction. 

Online supplemental material 

Fig. SI provides insight into the image analysis process (segmentation, 
threshold of detection, and newly visualized strains). Fig. S2 shows the 
comparison of our intensity measurements to published datasets. Fig. S3 
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shows abundance change events for each stress and provides the statisti- 
cal analysis for detection of significance. Fig. S4 shows detection of po- 
tential posttranscriptional regulation events by comparison to published 
microarray datasets. Fig. S5 provides examples for proteins that change 
abundance or localization during growth in H2O2 or nitrogen starvation 
and are essential during these growth conditions. Table SI summarizes 
intensity measurements and localization assignments under all conditions. 
Table S2 summarizes new assignments for protein localization. Table S3 
summarizes significantly changed proteins during starvation. Table S4 
summarizes significantly changed proteins in DTT. Table S5 summarizes 
significantly changed proteins in H2O2. Table S6 summarizes bimodal- 
ity distribution under stress. Table S7 summarizes localization changes 
under stress. Table 38 summarizes concerted movements of physically 
interacting proteins under stress. Table S9 summarizes newly acquired 
interactions of Accl in DTT. Table SIO summarizes yeast strains used in 
this study. Table 31 1 summarizes plasmids used in this study. Table 312 
summarizes primers used in this study. Table 313 summarizes systemati- 
cally removed proteins. Video 1 shows a fraction of the population of 
cells expressing mitochondrial protein Aim37-GFP in which it changes its 
mitochondrial localization to vacuole membrane under H2O2. Video 2 
shows a fraction of the population of cells expressing mitochondrial pro- 
tein Aim37-GFP in which it did not change its localization under H2O2. 
Online supplemental material is available at http://www.icb.org/cgi/ 
content/full/icb.201 301 1 20/DCl . Additional data ore available in the 
JCB DataViewer at http://dx.doi.org/10.1083/jcb.201301 120.dv. 
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